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Scaling of the conductances and the finite-size localization lengths is generalized to anisotropic sys- 
tems and tested in two dimensional systems. Scaling functions of isotropic systems are recovered 
once the dimension of the system in each direction is chosen to be proportional to the localiza- 
tion length. It is also shown that the geometric mean of the localization lengths is a function of 
the geometric mean of the conductivities. The ratio of the localization lengths is proportional to 
the square root of the ratio of the conductivities, which in turn is proportional to the anisotropy 
strength t, in the weak scattering limit. 
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Scaling theories have been successfully applied to the 
problem of Anderson localization |^,|^ involving the ef- 
fects of disorder on the nature of the electronic wave func- 
tion. The most remarkable result of the scaling theory is 
its prediction [|| that a continuous metal-insulator tran- 
sition exists in three dimensions, and all the states are 
localized in two dimensions, in the absence of magnetic 
field and spin-orbit interaction. The essential hypothe- 
sis of the one parameter scaling theory of localization is 
that the rate of change of the conductance when the size 
of the system changes is controlled by the conductance 
alone. The critical conductance Gc that separates true 
metals from insulators is estimated Q to be Gc = 0.1^. 
The scaling function should also be universal within a 
class that is characterized by a few general symmetries 
of the governing Hamiltonian. The scaling theory results 
are supported by a large number of numerical studies |^] 
in d=2 and d=3. Most notably, finite size scaling calcu- 
lations Q on the transmission properties of a quasi-one 
dimensional system explicitly demonstrated the existence 
of a universal scaling function close to the critical regime. 

Most of the previous work involves isotropic sys- 
tems. Recently, the problem of Anderson localization 
in anisotropic systems has attracted considerable atten- 
tion [||-|8|, largely due to the fact that a large variety of 
materials are highly anisotropic. It was recently shown 
that in a highly anisotropic system of weakly cou- 
pled planes, states are localized in the direction paral- 
lel and perpendicular to the plane at exactly the same 
amount of critical disorder, in support of the one parame- 
ter scaling theory which excludes the possibility of having 
a wave function localized in one direction and extended 
in the other two. However, several issues regarding the 
relation between the conductances in different directions 
were raised. Most importantly, the question of scaling of 
conductances and localization lengths was not resolved 

. Although anisotropy is known not to change the uni- 
versality and thus the critical behavior of the system [p| , 



the exact form of the scaling function, on the other hand 
is expected to depend on the anisotropy in the form of 
anisotropic physical parameters such as anisotropic hop- 
ping integrals or geometrical aspect ratios . 

Extending the scaling argument to an anisotropic sys- 
tem, we assume that the logarithmic derivative. Pi, of the 
dimensionless conductance gi, in any direction will be a 
function of the conductance in that direction as well as 
other directions, 



dloga 



m9^}) 



(1) 



where a is an appropriate length scale. All the gi become 
relevant scaling parameters. All other physical quanti- 
ties, such as anisotropic hopping integrals or anisotropic 
geometrical shapes, should enter only through the con- 
ductances gi. Exactly the same argument can be applied 
to the scaling function of localization length, obtained 
from transfer matrix calculations with a quasi-one di- 
mensional geometry of cross section Mj x Mk, 
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where Xi is the finite size localization length in the di- 
rection i and £,i {I = 1,2,3) is the localization length for 
the infinite system. The fundamental assumption in Eq. 
(2) is that localization lengths provide the only charac- 
teristic length scale. Once the characteristic lengths are 
measured in terms of the localization lengths in the cor- 
responding directions, the scaling behavior of the system 
within the same universality class, are governed by the 
same equation. 

The scaling functions Pi and h describe the behavior of 
both systems with isotropic Hamiltonians but non-cubic 
geometry as well as systems with anisotropic Hamilto- 
nians. Scaling in anisotropic systems in general is not 
known. Only when the conductances in all directions are 
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the same, then the scaUng function (Si will be exactly 
the same as that of a cubic isotropic system. For an 
anisotropic system, this can only be achieved by choos- 
ing an appropriate geometry which might not be known 
a priori. As an example we will see that indeed such a 
procedure works in a system with highly anisotropic hop- 
ping. We will demonstrate that in two dimensional sys- 
tems, Eq.(2) can be applied straightforwardly such that 
all the data are described by the scaling functions of the 
isotropic system. Furthermore, we will also show that the 
geometric mean of the localization lengths is a universal 
function of the geometric mean of the bare conductivity, 
and their ratio can also be estimated in the weak scat- 
tering limit. These results follow directly from applying 
the basic idea of scaling theory. 

We consider the following Hamiltonian for an 
anisotropic 2d disordered model 
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where n labels the sites of a square lattice. The on-site en- 
ergies en are independently distributed at random, within 
an interval of width W. The second term is taken over all 
pairs of nearest neighbor (n.n.) sites, and the hopping in- 
tegral i„m=l or t (<1), depending on hopping directions. 
As a convention, we have assigned the direction with the 
large {tnm = 1) and the small (t„m = t) hopping value 
as the parallel (||) and the perpendicular (_L) directions, 
respectively. 

In two dimensional systems, Eq. (2) can be written as 
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where f(x)=h(x)/x is the scaling function for isotropic 
systems. We have used the transfer matrix method ||] 
to calculate the finite size localization length Xi{Mj) for 
many Mj = 1,2) (M=24, 48, 96, 120, 150, 300) and 
W=2-14 and several t and E, for both directions. Figure 
1 shows that all of our raw numerical data for both aJ[,j 
and A^^ for different anisotropics t different disorder W 
and different energies E follow one universal curve, by 
appropriately choosing the localization length in the two 
directions, ^|| and The solid line through the data 
in Fig. 1, is the 2d isotropic scaling function. This is a 
direct confirmation of the scaling relation Eq. (4) . 

An important consequence of Eq. (4) is that at the 
critical point, if any, the geometric mean of the ratio 
of the finite-size localization length to the cross-section 
width is a constant. This was indeed found to be true 
but interpreted instead as a result of possible conformal 
invariance. We point out that at the critical point, the 
geometric mean of the conductances along the different 
directions may not be a constant. This behavior of the 
conductances is different from that of \m/M and needs 
further study for its complete understanding. 

To further test the scaling idea we have calculated the 
conductance G in the two different directions for our 



anisotropic system. From the Multichannel Landauer 
formula ||ll|,@, G = ^Tr{th), where t is the trans- 
mission matrix. With anisotropic hoppings, one should 
choose a geometry other than the square such that the 
conductance is the same in all the directions and then 
scale up the size of the system The conductance 
should remain isotropic if one parameter scaling theory 
is correct We have tested this idea in a 2d system 
with t = 0.1. The ratio of the two localization lengths 
was found to be 10 at W = 3.6. We have scaled up the 
system of a rectangle of size Mx N by a factor 4, and 
from Fig. 2, one can clearly see that although the con- 
ductance becomes extremely small it remains isotropic, 
in agreement with the predictions of the one parameter 
scaling theory For a square geometry and with the 
same parameters as in Fig. 2 , the conductances in the 
two directions would diverge rapidly as the system size 
scales up. 

Another length rescaling aspect can be seen by con- 
sidering the self-consistent theory of localization. It was 
shown in an earlier work that in order the localiza- 
tion criteria to be the same in all directions, the length 
scale has to be chosen proportional to the square root of 
the bare conductivity. This leads to an equation for the 
metal-insulator transition that is exactly the same as that 
of the isotropic system except that both the bare conduc- 
tivity and the effective lattice constant are replaced by 
their geometric means. A direct consequence of this for- 
mulation is that the geometric mean of the localization 
(or correlation) lengths should be only a function of the 
geometric mean of the bare conductivities, ie, 
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where <>g denote the geometric mean of the values in 
the two directions, do is the bare conductivity and Sf is 
the Fermi surface area that enters through the relation 
(To ~ Sf£. i is the mean free path. // is a function that 
can be obtained via the Potential Well Analogy(PWA) 
or the self consistent theory of localization. Using the 
PWA Ref.l3 obtained ^ = 2.72iexp[n^hao/e^]. 

Eq. (5) can be easily checked in the weak disorder 
limit, at which the geometric mean of the bare conduc- 
tivity can be shown to be (Tq = IdV^t / ttW^ , within the 
Coherent Potential Approximation (CPA) [14]. For the 
2d anisotropic system, Sf{E = 0) = 47r\/TTt^. We 
have plotted Sf < >^/^ versus Iby^t/nW^ , and 

find that the data fall into one universal curve for all 
the different anisotropics t and disorder W. This weak 
scattering limit behavior of the geometric mean of the 
localization lengths versus the geometric mean of the con- 
ductivities is very suggestive of the way the localization 
lengths have to scale. The full expression, valid for all 
disorder strength. 
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can also be evaluated. S = Si — iTj2 is the self-en ergy 
obtained by solving a self-consistent equation |]5|, |l3| , |l5| . 
This shows remarkable good scaling as shown in Fig. 3, 
including results for E=0, as well as for E 0. The 
curve in Fig. 3 shows how the geometric mean of the 
localization length depends on the geometric mean of the 
bare conductivity in a universal fashion, independent of 
the anisotropy, energy and disorder. These results are 
a strong confirmation of scaling in anisotropic systems. 
Notice that the geometric mean of the conductivities < 
(To >g, and not of the conductances, is the appropriate 
quantity that gives the same results as in the isotropic 
case. It is therefore appropriate that < uo >g will be used 
in the interpretation of experiments in highly anisotropic 
systems. 

The ratio of the localization lengths can be obtained 
by carrying the length rescaling idea further. We can see 
that the conductances in all the directions should be the 
same, if the dimension of the system is proportional to 
the localization length in that direction. This implies the 
following relation 

ai is the exact value of the conductivity, aio is the bare 
conductivity which can be calculated within the CPA, ai 
is the correction of the bare conductivity in the i direc- 
tion. It is very difficult to calculate the correction ai, but 
it approaches one in the weak scattering limit. In Fig. 
4, we show the results of versus 1/W for different 

anisotropics t and energies E. In the weak disorder limit, 
we can approximate cr by CTq, and this is shown as open 
symbols in Fig. 4. Notice that in the weak disorder limit, 
W 0,and for t ^ 0, (Jo±/'^o\\ ~ and by using Eq. 
(7), one obtains for the ratio of the localization lengths 
, i±/i\\ ~ t- This behaviour is clearly seen in Fig. 4 
for large 1/W. Agreement with the CPA results for the 
conductivity are excellent for weak disorder. Deviation 
of the ratio from the open symbols for strong disorder 
indicates that the true conductivity at length scale ^ is 
stongly normalized compared with the bare conductivity. 
However, it is notable that the trend of the ratio as W 
increases is captured by the simple expression. For large 
W, no dependence on E should be expected for small E, 
thus the ratios converge to the same value for different E 
with t=0.3, as can be seen in the insert in Fig. 4. 

In summary, we have performed an extensive numeri- 
cal study of the scaling properties of highly anisotropic 
systems. Scaling functions of isotropic systems are recov- 
ered once the dimension of the system in each direction 
is chosen to be proportional to the localization length. In 
the localized regime, the ratio of the localization lengths 
is proportional to the square root of the ratio of the con- 
ductivities which in turn is proportional to the strength 
of the anisotropy t (i.e. £,±/£,\\ ~ t). Recall that in the ex- 
tended regime the ratio of the correlation lengths is 
proportional to the ratio of the conductivities (i.e. 



= coi/cToll 1/^^ )• It was also shown that the geo- 
metric mean of the localization lengths is a function of 
the geometric mean of the conductivities, not of the con- 
ductances. Finally, it was numerically shown that the 
conductances along the two different directions of the 
anisotropic system are the same, provided that the di- 
mension of the anisotropic system is proportional to the 
localization length in this direction. This procedure can 
be easily used in other anisotropic systems. 
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FIG. 1. The numerically determined scaling function for 

the 2D anisotropic system for different anisotropic constants 
t, different energies E and disorder W. The solid lino through 
the data is the 2D isotropic scaling function. The y-£ixis is 
^jXi{Mj)/^iMj, while the x-axis is Mj/^j. The index i and 
j can be either the parallel or the perpendicular direction, 
respectively. 

FIG. 2. The conductance G in units of e'^/h of an 
anisotropic system M x N, versus M for t=0.1 and E=0. No- 
tice that G along the two directions is exactly the same. 

FIG. 3. The product of the Fermi surface Sf with the 
geometric mean of the localization lengths < 5±C|| >^^^ is 
plotted versus the geometric mean of the bare conductivities 
< Cox foil >^/^ for all the energies E, t's and W. 

FIG. 4. The ratio of Cx/^^n is plotted versus 1/W, for 
t=0.1, 0.3 and 0.6 with E=0. and are the localiza- 
tion lengths along the two propagating directions. The solid 
symbols arc the numerical results, while the open symbols are 
the CPA results. In the insert the numerical results of ^±/t^\\ 
versus 1/W is plotted for t=0.3 with E=0.0, 1.5 and 2.0. 
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